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The vorticity development is studied in the reaction plane of peripheral relativistic heavy ion 
reactions where the initial state has substantial angular momentum. The earlier predicted rotation 
effect and Kelvin Helmholtz Instability, lead to significant initial vorticity and circulation. In low 
viscosity QGP this vorticity remains still significant at the time of freeze out of the system, even if 
damping due to the explosive expansion and the dissipation decreases the vorticity and circulation. 
In the reaction plane the vorticity arises from the initial angular momentum, and it is stronger than 
in the transverse plane where vorticity is caused by random fluctuations only. 

PACS numbers: 89.65.-s, 89.75. Da, 05.45.Tp. 



I. INTRODUCTION 

The strongly interacting quark gluon plasma (QGP) [T] 
has raised many interesting questions about the physics 
of ultra dense hot matter produced in relativistic heavy 
ion collisions. Many transport and hydrodynamic mod- 
els [5] are popular in analyzing the evolvement of the 
viscous QGP and its properties. In peripheral heavy ion 
reactions due to the initial angular momentum the ini- 
tial state of the fluid dynamical stage of the collision dy- 
namics has a shear flow characteristics, and this leads to 
rotation [3] and even Kelvin Helmholtz Instability (KHI) 
[4] in the reaction plane, for low viscosity Quark-gluon 
plasma. This possibility was indicated by high resolu- 
tion Computational Fluid Dynamics (CFD) calculations 
using the PIC method. We study the development of 
these processes in 3+1 dimensional (3+1D) configuration 
to describe the energy and momentum balance realisti- 
cally. 




FIG. 1. The sketch of a collision. Figure (a) is in the 
transverse, [x-y], plane and (b) is in the reaction, [x-z] plane. 
The almond shape in the middle of figure (a) is the participant 
zone of the event. Right after the collision, streaks are formed 
and the top streaks move along the z direction while bottom 
ones move along the — z direction. This velocity shear may 
lead to the Kelvin Helmholtz Instability, a wave formation, 
on the interface between the top and bottom sheets. 



In idealized 2+1D model calculations the dissipation 
due to 3D expansion is neglected, and thus the dissipation 



due to the 3D viscous expansion is also neglected, which 
results in unrealistic estimates. The vorticity of the flow 
is especially sensitive to such over-simplifications: while 
in 2+1D the integrated vorticity is conserved in perfect 
fluid flow, the decrease of vorticity is essential in realistic 
3+ ID CFD description. 




FIG. 2. (Color online) The 3-dimensional view of the col- 
lision a few fm/c after the impact corresponding to the situ- 
ation illustrated in Fig. [l](a). The projectile spectators are 
going along the z-direction, the target spectators are going 
along the —z axis. We assume that the participants in the 
middle, form a cylinder with an almond shape profile and 
tilted end surfaces, where the top side is moving to the right 
and the bottom is moving to the left. The participant cylinder 
can be divided into streaks, each streak has its own velocity, 
as shown in Fig. [I] (b) . The velocity differences among the 
streaks result in the KHI effect. 

In ref. [5] the angular momentum is assumed to have 
significant effects on the longitudinal flow velocity and on 
its distribution in the transverse plane, so that it gives 
rise to vorticity and polarization. The arising polariza- 
tion is also studied in [BJ, where a laminar shear flow is 
assumed where each layer has different velocity, which is 
quite similar to the initial state velocity profile depicted 
in Fig. [T] This type of initial state is described in great 
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detail in [7]. In our present fluid dynamical calculation 
we use this initial state model, which is tested in several 
model calculations in the last decade. It describes cor- 
rectly the initial shear flow characteristics. The initial 
angular momentum is based on the assumption that the 
initial angular momentum of the participants (based on 
straight propagation geometry) is streak by streak con- 
served. This leads to strong shear flow, shown in Fig. 
[2j This model assumes that the incoming Lorentz con- 
tracted nuclei interpenetrate each other, and the lead- 
ing charges in each streak are slowed down by the large 
string-rope tension. This takes about 3-5 fm/c for heavy 
nuclei depending on the impact parameter. Then local 
equilibration is reached and the fluid dynamical evolution 
starts. 

In this work we study the development of vorticity in 
high energy heavy ion reactions, and the development of 
the above mentioned specific effects, which may arise in 
low viscosity [8, 9 QGP. 

We also compare the classical vorticity characteristics, 
and the effect of the most dominant relativistic general- 
izations. These may provide an insight into the possi- 
bilities of using these effects to precision studies of the 
transport properties of QGP. 

It is also important to mention that the shear flow 
velocity profile is essential from the point of instability. 
According to our initial state description of the fluid dy- 
namical calculations (Fig. 4 of [1]) the velocity profile 
along the x-axis is not linear but has ai~ tan(w z ) shape 
(dotted line in Fig. [3]), which may lead to KHI, in con- 
trast to a x ~ arctan(w z ) shape, which does not. 
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FIG. 3. The velocity along the z axis, v z , represented by the 
dotted curve is calculated in our CFD model and presented in 
Fig. 4 of Ref. [3]. Here we have to mention that the velocity of 
the dotted curve will induce the KHI effect, while the dashed 
curve will not, see chapter 8 of Ref. [TU] . 



In the following we will discuss recent approaches to 
flow vorticity in high energy heavy ion collisions, and 
present the vorticity development and its distribution in 
the reaction plane. 



In classical physics for incompressible, perfect fluids 
vorticity exhibits an impressive conservation law: the 
conservation of circulation. In high energy heavy ion 
physics the vorticity definition must be modified, and 
the mentioned idealizations are not applicable to ener- 
getic heavy ion reactions. Still, as CFD calculations in- 
dicate typical flow patterns and instabilities may occur 
here also. Thus their studies can provide insight into the 
properties of the QGP fluid. 
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FIG. 4. (Color online) The classical (a) and relativistic 
(b) weighted vorticity, Q zx in units of c/fm, calculated in the 
reaction, [x-z] plane at t=0.17 fm/c after the start of fluid dy- 
namical evolution. The collision energy is ^/sjvjv = 2.76 TeV 
and b = 0.7b max , the cell size is dx — dy = dz = 0.4375 fm. 
The average vorticity in the reaction plane is 0.1434 / 0.1185 
c/fm for the classical / relativistic weighted vorticities respec- 
tively. 



II. THE VORTICITY 

a. Classically, in the reaction plane, [x-z] , the vor- 
ticity is defined as: 

uj y = w xz = -lo zx = -(d z v x -d x v z ) (1) 

where the x, y, z components of the 3-velocity v are de- 
noted by v x , v y , v z respectively. In this definition we 
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have already included the factor \ for the symmetriza- 
tion to have the same magnitude of vorticity as for sym- 
metrized volume divergence or expansion rate. 

Here we study the vorticity in the reaction plane (in the 
[x-z]-plane), at different times in the CFD development. 

In 3-dimensional space the vorticity can be defined as 

u = -rot d = -Vxd (2) 
2 2 w 

(but this cannot be generalized to four or more dimen- 
sions). The circulation of the flow is the integral of the 
velocity along a closed curve, C, with the line element 
dl. It is defined as 

T = * v dl = j 2ui dA , 

where A is an (arbitrary) surface surrounded by the curve 
C, and the normal of its surface element is dA. For 
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FIG. 5. (Color online) The classical (a) and relativis- 
tic (b) weighted vorticity, Q, zx , (c/fm), calculated in the 
reaction [x-z] plane at t=3.56 fm/c. The collision energy 
is ^/sjvjv = 2.76 TeV and b = Q.7b max , the cell size is 
dx — dy = dz — 0.4375/m. The average vorticity in the 
reaction plane is 0.04845 / 0.07937 c/fm for the classical / 
relativistic weighted vorticity respectively. 

rotationless (potential) flow the circulation vanishes. In 
peripheral heavy ion reactions the flow is rotational. In 
classical fluid dynamics the circulation is constant in a 



fluid along the line of motion of a fluid element if the 
viscosity of the flow vanishes and flow is barotropic, i.e. 
the pressure depends only on the density of the fluid, 
P = P(p). These conditions are not satisfied for QGP in 
heavy ion reactions. 

In the relativistic case the formulation of fluid dynam- 
ics is more involved, and the formalism must be modified 
[TT1 H2] . The three- velocities are replaced by four- veloc- 
ities, the derivatives should take into account the world 
lines of the particles and the changes of these with time 
and the mass density is replaced by the energy density 
where the pressure will have a non-negligible role. This 
last modification e.g. leads to a modified relativistic defi- 
nition of the vorticity and circulation, which includes the 
specific enthalpy of the fluid [13] ■ This would then ex- 
tend the validity of the conservation of circulation under 
the same conditions. On the other hand this modified 
relativistic vorticity would then have a different dimen- 
sion and it would not be conserved anyway in heavy ion 
reactions where the conditions of validity are not satis- 
fied. 
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FIG. 6. (Color online) The classical (a) and relativis- 
tic (b) weighted vorticity tt zx , (c/fm), calculated in the re- 
action [x-z] plane at t=6.94 fm/c. The collision energy 
is ^/snn = 2.76 TeV and b = 0.7b max , the cell size is 
dx — dy — dz — 0.4375/m The average vorticity in the 
reaction plane is 0.01555 / 0.05881 c/fm for the classical / 
relativistic weighted vorticity respectively. 

Thus, in the present work we intend to keep the possi- 
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bility of comparison to the classical vorticity and circu- 
lation, in a relativistic system where the (baryon charge) 
density and temperature change violently, the pressure 
is not negligible and does not depend on the density 
only. Still the energy and momentum, (T°° and T lQ ) 
are strictly conserved, the total energy remains constant 
while the momentum vanishes (in the CM. frame). 

Keeping in view fluid dynamics with fluid elements, 
we introduce a weight proportional to the energy con- 
tent of the fluid element. The energy distribution (and 
thus the local angular momentum) of the matter is highly 
non-homogeneous in a heavy ion reaction. To reflect our 
physical situation better, we weight the contribution of 
our fluid cells by the local energy density in the reaction 
plane, however, without changing the average vorticity 
in the layer of the reaction plane. Thus we weight with 
a distribution normalized to unity. We define an energy 
density weighted, average vorticity as 



Q zx = w(z,x) 



(3) 



so that this weighting does not change the average cir- 
culation of the layer, i.e., the sum of the average of the 
weights over all fluid cells is unity, (w(z, x)) = 1, both 
if we consider one y-layer only, or if we consider all y- 
layers. This weighting does not change the average vor- 
ticity value of the set, just the cells will have larger weight 
with more energy content, T 00 . 

As we have discretized fluid cells, we study separately 
the reaction plane, one y-layer, or the whole system, all 
y-layers. The total energy content of a cell at point (z, x) 
or the corresponding i,k, is Eik = T Q0 (z,x). The total 
energy in a y-layer (or in all y-layers) is E tot = J2 ik E ik, 
while the number of the cells is in a y-layer (or in all 
y-layers) is N ce u . Thus the average energy for for a fluid 
cell is E to t I ^ ceil in both cases. We divide the actual en- 
ergy of a fluid cell, E ik , with the average fluid cell energy 
Etot I ^ ceil: so that the vorticity values on average will 
remain comparable with the non-weighted values, but 
still larger energy cells will have more weight. The to- 
tal energy ^2 ah ceils T®® remains exactly constant in our 
case (in the numerical calculation to 10~ 16 accuracy). Of 
course this is not true for a single layer like the reaction 
plane. 

Our weight then, should be proportional with the local 
energy density 



E. 



ik 



(E to t/N ceU ) 



Within the reaction plane, at a given moment of time, 
the cells carrying larger amount of energy will get a larger- 
weight than those, which carry less energy. The edge 
cells, carrying less energy, show stronger fluctuations. 

This average vorticity weighted by the cell energy, 
which is a conserved quantity for all cells, provides a 
possibility to compare the results to classical systems and 
their features. 




FIG. 7. (Color online) The classical (a) and relativistic 
(b) weighted vorticity calculated for all [x-z] layers at t=0.17 
fm/c. The collision energy is ^/sjvjv = 2.76 TeV and b = 
0.7b m ax, the cell size is dx = dy = dz — 0.4375 fm. The 
average vorticity in the reaction plane is 0.1971 / 0.19004 c/fm 
for the classical / relativistic weighted vorticity respectively. 



An alternative method to present the relativistic vor- 
ticity is presented in ref. [T3] by weighting with the spe- 
cific entropy, would provide conserved circulation, T, if 
the pressure were exclusively density dependent. As it 
was mentioned earlier, this is not the case for quark-gluon 
plasma, so for us this advantage is not realized, while this 
weighting changes both the value and the dimension of 
the vorticity, so it would make the comparison to classical 
results difficult. 

The weighted vorticity in the reaction, [x-z] plane at 
different time steps are shown in the Figs. [4j [5j [6] 

b. In relativistic flow, we are following the defi- 
nition of in Ref. [T2], for the relativistic case. The ex- 
pansion rate, O and the vorticity tensor, w£f are defined 
as 



e 



(4) 



(5) 



were for any four vector g M the quantity V a <? p 



A£ dpq» = AP (ft and A" 



fr 



u^u v . This leads 
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FIG. 8. (Color online) The classical (a) and relativistic 
(b) weighted vorticity calculated for all [x-z] layers at t=3.56 
fm/c. The collision energy is ^/sjvjv = 2.76 TeV and b = 
0.7 bmax, the cell size is dx = dy = dz — 0.4375 fm. The 
average vorticity in the reaction plane is 0.0538 / 0.10685 c/fm 
for the classical / relativistic weighted vorticity respectively. 



FIG. 9. (Color online) The classical (a) and relativistic 
(b) weighted vorticity calculated for all [x-z] layers at t=6.94 
fm/c. The collision energy is ^snn ~ 2.76 TeV and b = 
0.7b ma;c , the cell size is dx = dy = dz — 0.4375 fm. The 
average vorticity in the reaction plane is 0.0159 / 0.05881 c/fm 
for the classical / relativistic weighted vorticity respectively. 



to 

= 2 [( d » u " - d " u ») + (u"u a d a u v - u u u a d a u»)] 
= X - [(cW - d"u u ) + (u»d T u u - u„d T u»)] , 

where drU^ = u M = v^daU 11 is the proper time derivative 
of u". 

In our initial state in the middle part of the collision 
the streaks are stopped after two Lorentz contracted ions 
collide and interpenatrate each other. If the acceleration 
of the fluid elements is negligible compared to the ro- 
tation, |<9 r it M | -c |9 K it z |, this assumption holds for the 
initial moments in our model, which has strong initial 
shear flow. Thus the second term can be dropped. This 
is also the case considered in Ref.|llj. while this work is 
studying the vorticity in the [x-y] plane instead of the 
reaction plane. In this case the relativistic vorticity is: 

= ~ d n u v) ( 6 ) 



where 

d v = (d ,d x ,d y ,d z ) 

Ufj, = 7(1, -V X , ~Vy, ~Vg) 

Let us expand Eq. (16]) in four dimensions: 
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and we see that the vorticity is an antisymmetric tensor. 
Here for the vorticity development in the reaction plane 



G 



we calculate LO24' 



= ^(d z -fv x - d x jv z ) 

= ]p(d z v x - d x v z ) + ^(v x d z ~/ - v z d x j) (8) 

The fluid cells were weighted the same way as in the 
non-relativistic vorticity estimate. Due to the fact that 
the relativistic vorticity includes the relativistic 7 factor, 
eq. Q, the vorticity increases, especially at the edges 
where the cells have larger flow velocities. 

The weighted relativistic vorticity distributions in the 
reaction, [x-z] plane at different time steps are shown 
in the Figs. |4j [5j [6] At the last time step presented, 
in the reaction plane we have already an extended area 
occupied with matter. In case of peripheral reactions the 
multiplicity is already small, thus the fluctuations in the 
reaction plane are considerable. In the relativistic case 
the outside edges show larger vorticity and the random 
fluctuations are still strong. 

III. WEIGHTED VORTICITY DISTRIBUTIONS 

The relativistic vorticity distributions have increased 
amplitudes compared to the classic ones, due to the rela- 
tivistic 7-factor and its derivatives in the relativistic ex- 
pression. This is especially visible at the edges where the 
flow velocities are the largest. 

The amplitude of weighted vorticity decreases with 
time as the matter expands. Random fluctuations are 
apparent at late times for the dilute matter, especially at 
the edges, and particularly for the relativistic vorticity, 
which has enhanced amplitudes. 

With increasing time the fluid expands, and outside 
more dense shell develops with a less dense central zone. 
This feature is also apparent in the vorticity distribution 
of the mater in the central, reaction plane, where the 
central part has smaller weighted vorticity amplitude. 

The non-central, parallel layers, at increasing y-values, 
have similar positive weighted vortices, VL y zx . These out- 
side layers have narrower boundaries in the ^-direction 
and these fall into the denser outside zone of the expand- 
ing matter. 

We can average the weighted vorticity distributions 
over all layers parallel to the reaction plane. This com- 
pensates the low density central zone in the reaction 
plane and leads to a more uniform, layer-averaged distri- 
bution, with higher positive peak amplitudes and smaller 
negative values. The last presented time-step at around 
t = 7 fiji/c is still strongly fluctuating for the relativis- 
tic case. The weighted vorticity averaged over all layers 
parallel to the reaction, [x-z] plane at different time steps 
are shown in Figs. [7J |8j [9j The dominant effect of the 
relativistic treatment is visible the most in Fig. [9] Com- 
paring eqs. [2] and [8] the role of the relativistic 7 factor is 
apparent. As discussed after eq. ^ the large peripheral 



velocities make the vorticity large at the external surfaces 
of the expanding system. At these external regions the 
matter density is small and the relative density fluctua- 
tions are large, so the fluctuations of the surface region 
vorticities are large. 

IV. THE PROPERTIES OF VORTICITY AND 
CIRCULATION 

In section [TT] we introduced the classical vorticity and 
circulation with the conservation laws for the circulation 
in case of certain conditions. These can be extended to 
the relativistic case if we define the relativistic circulation 
(chapter 14.3) as: 

rrc) = — I wu^Sxn, 

m Jc 

where the weight, w, is the specific enthalpy, e + P, over 
the baryon number, n, and m is the effective mass per net 
baryon. Applying the Euler equation for perfect fluids 
results in 

dT(C) 1 I d v P If 

» = — 9 5x v + — <p wu^Sun , (9) 

Ot m Jc n m J c 

where the last term vanishes along the flow stream line 
because u^du^ — 0. Thus, if the pressure, P depends 
on the density, n, only the second term for a closed loop 
integral also disappears, and the circulation remains con- 
stant just as in the classical case. 

In our problem these conditions are not satisfied, so 
we used the same weighting as in the classical case, and 
evaluated the circulation with the weighting as in the 
classical case to enable us the comparison, see Fig. [T0| 
We performed the integral over the [x-z] surface for the 
weighted relativistic vorticities avaraged over all y-layers. 

Let us take the surface area of the reaction plane, which 
is filled with fluid and take the bounding curve at the out- 
side edges of the fluid in this plane. This curve expands 
with the fluid. In our case in rapidly expanding and non- 
perfect fluid the circulation as well as the vorticity in in 
the [x-z] reaction plane are both decreasing. 

For us the rate of this decrease is important, to see if 
we can still detect the vorticity and circulation at freeze 
out. Notice that we calculated only the Q zx component 
of the weighted vorticity distributions. Due to the close 
to spherical expansion, the direction of vorticity may de- 
velop into different directions. This also contributes to 
the decrease of the circulation. 

We calculated and presented the weighted vorticity dis- 
tribution in the reaction plane, and then also calculated 
it for all the [x-z] layers at different y-values, and took the 
average of these vorticity distributions, see Figs. [7J [5J |9j 
[TT]3, [l2p. The overall vorticity is positive, this originates 
from the initial shear flow configuration. It decreases 
with time, not only because the weight of the reaction 
plane decreases due to the expansion, but also because 
of the viscous dissipation in the 3D expansion. 
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FIG. 10. The time dependence of classical circulation, F(t), in units of [fmc], calculated for all [x-z] layers and then taking 
the average of the circulations for all layers. The collision energy is ^/snn = 2.76 TeV and b = Q.7b max , the cell size is 
dx — dy = dz = 0.4375 fm (left). For comparison another initial state configuration was also tested for the same collision 
energy but b = 0.5 b max , the cell size is dx — dy — dz = 0.585 fm (right). This configuration shows also the rotation, but 
due to its less favorable parameters it does not show the KHI. Although at this impact parameter, which is less peripheral the 
reaction plane has a larger area filled with matter, nevertheless the initial classical circulation is less by about 15%. For the 
more peripheral case with smaller numerical viscosity the circulation decreases with time faster and the circulation for the two 
cases becomes equal around t = 8 — 10 fm/c. 



The examples above are for the configuration where 
the Kelvin-Hclmholtz Instability (KHI) is predicted to 
occur in heavy ion collisions at LHC. 

We also studied, how the vorticity and circulation 
changes with increased (numerical) viscosity and for more 
central collisions, where the occurrence of KHI is not pre- 
dicted by the CFD model, see Figs. [U^,[l2^. 

In ref. [1] it was suggested that the KHI leads to an 
increased v\{y) peak. We search for other, preferably 
more sensitive and more specific experimental methods 
to identify Rotation and KHI, and possibly also separate 
the two effects. 

The vorticity is also strikingly different for the configu- 
rations which are adequate for KHI, and for those which 
not, and show only rotation. The initial and intermediate 
time stages are compared for two different configurations 
where KHI may and may not appear. 

The primary reason of the difference with KHI is that 
in the highly peripheral reactions the profile height is 
smaller but the asymmetry between the target and pro- 
jectile contributions is larger at the edges, and thus the 
shear in the matter is considerably larger. Consequently 
both the maximum value of vorticity and the average are 
larger for the more peripheral configuration: the initial 
average vorticity is almost 3-times larger in the more pe- 
ripheral configuration favourable for KHI than in the less 
favorable one. 

At the later time, t = 3.56 fm/c, the difference is still 
large, the KHI formation leads to a vorticity, which is 



about a factor two larger, Fig. 12 



This, strong dependence on impact parameter arising 
in low viscosity matter is very promising from the point 
of view of the observability of the effect. 

The possibility for observation mentioned in ref. [4], 
by the position of the rotated v\{y) peak, is also visible 
in Fig. [12] where the vorticity peak position is seen in the 
two configurations. This peak position may be coupled 
to the peak position of the earlier mentioned wi-peak. 
Here in the configuration favorable for KHI the forward 
rotation angle of the peak is 37°, while in the less favor- 
able configuration it is 33°. Thus, the KHI arising only 
in low viscosity matter leads to a special increase in the 
rotation and vorticity. The observability via finding the 
motion of the collective (y-odd) v\ (y) peak is not easy in 
present LHC experiemnts because the initial state fluctu- 
ations contribute to large longitudinal cm. fluctuations, 
apart of azimuthal fluctuations in the transverse plane 



V. CONCLUSIONS 

An analysis of the vorticity and circulation develop- 
ment was performed for peripheral Pb+Pb reactions at 
the CERN LHC energy of v / i]^ = 2.76 TeV. The initial 
peak vorticity was more than 10 times larger (exceeding 
3 c/fm) than the one obtained from random fluctuations 
in the transverse plane, of about 0.2 c/fm [11]. The rea- 
son is in the high initial angular momentum arising from 
the beam energy in non-central collisions. 




z(fm) 

FIG. 11. (Color online) The weighted relativistic vorticity 
calculated for all [x-z] layers at t=0.17 fm/c. The collision 
energy is ^/sivjv = 2.76 TeV, b = 0.5 fo ma a; for (a) and b = 
0-7b m ax for (b). The configuration (a) is not favoring KHI 
while (b) is. The cell size is dx = dy = dz = 0.585/0.4375 
fm. The average weighted vorticity in the reaction plane is 
0.07241/0.19004 c/fm for the two cases respectively. Notice 
the different color coding scales for (a) and (b). 



Although the vorticity and circulation decreases 
rapidly due to the explosive expansion of the system still 
at 4 fm/c after the beginning of fluid dynamical expan- 
sion the peak vorticity is above 3 c/fm in favorable con- 
figurations with KHI development and it reaches up to 1 
c/fm at the same time for less favorable configurations. 

This makes it promising to observe the consequences 
of this rotation and its sensitivity to turbulent configu- 
rations. Not only the peak vorticity exceeds earlier esti- 
mates from random fluctuations, but also the predicted 
average vorticity is substantial, it reaches 0.2 c/fm at 
favourable initial configurations with KHI, while the av- 
erage vorticity originating from random fluctuations is 
vanishing [TTJ [T3j . The estimated angular deflection aris- 
ing from this random origin "chiral vortaic effect (CVE)" 
[13], is small, cos(<5</>) ~ 10~ 5 . The effect arising from 
the initial angular momentum of a peripheral collision is 
bigger. This improves the observability of the average 
vorticity compared to CME and CVE predictions. 

In this work we just repeat the earlier mentioned ob- 
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FIG. 12. (Color online) The weighted relativistic vorticity 
calculated for all [x-z] layers at t=3.56 fm/c. The collision 
energy is ^Jsnn = 2.76 TeV, b = 0.5 b max for (a) and b = 
0.7b m ax for (b). The configuration (a) is not favoring KHI 
while (b) is. The cell size is dx = dy = dz = 0.585/0.4375 
fm. The average vorticity in the reaction plane is 0.05242 
/ 0.10685 c/fm for the two weighted vorticity respectively. 
Notice the different color coding scales. 



servable signatures related to the (rapidity odd compo- 
nent) of the directed, v\ flow, which is a promising pos- 
sibility for the observations [TS]. The time sequences of 
the results show that the maximum of vorticity in the 
side regions rotates, with time it moves forwards on the 
top, the maximum reaching the positive z-side for the lat- 
est times presented. The even flow harmonics are much 
less sensitive to this change. The rapidity width of the 
V2{y) may be effected weakly. Other methods in different 
correlation observables are more directly connected to ro- 
tation in the flow and will be addressed in forthcoming 
publications, (e.g. [d]). 

The question arises, how the surface energy influences 
the rotation and the KHI. The external surface of the ex- 
panding QGP is significant, as in the interior the quark 
gluon fluid has weak interaction (asymptotic freedom) 
and small viscosity. However, the surface energy has the 
strongest effect on the hadronization of QGP as first de- 
scribed in ref. [TJj. On the other hand the collective 
flow, as well as the KHI, develops primarily in the early 
QGP phase as the approximate quark number scaling in- 
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dicates. At RHIC and LHC energies, at this early stage 
the surface energy is negligible compared to the energy 
of the collective flow, and it primarily hinders the early 
emission of low energy hadrons from the plasma, while 
it does not hinder the rotation. In the KHI the situation 
is more involved. In the cases where the KHI develops 
between two fluids (e.g. water/air or oil/air) the large 
surface tension (e.g. between oil and air) can hinder the 
development of KHI, as well known to sailors for cen- 
turies. If the KHI develops within one fluid (like in air or 
QGP), due to the large shear between the two fluid lay- 
ers, there is no surface tension in the conventional sense, 



but the layer with the high shear may have extra energy, 
and so it can lead to an effective surface tension, which 
can weaken the development of KHI also. See estimates 
inref. [IS]. 
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